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Abstract. The one-dimensional half-filled Hubbard model is considered at zero 
temperature within the cellular dynamical mean-field theory (CDMFT). By the 
computation of the spectral gap and the energy density with various cluster and bath 
sizes we examine the accuracy of the CDMFT in a systematic way, which proves 
the accurate description of the one-dimensional systems by the CDMFT with small 
clusters. We also calculate the spectral weights in a full range of the momentum for 
various interaction strengths. The results do not only account for the spin-charge 
£H • separation, but they also reproduce all the features of the Bethe ansatz dispersions, 

Q | implying that the CDMFT provides an excellent description of the spectral properties 

i of low-dimensional interacting systems. 
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1. Introduction 



Strong correlations arising from the electron-electron interactions have been one of the 
key issues to solve in the condensed matter physics. The study on the strong correlations 
was initiated for the description of interacting electrons in localized orbitals. It has 
been accelerated by the discovery of interesting correlation properties such as high- 
temperature superconductivity or the colossal magnetic resistance. 

The main difficulty in the study of electron correlations is that even the simplified 
theoretical models, which include only a small number of important degrees of freedom, 
are very challenging to solve. The advent of the dynamical mean-field theory (DMFT) [1] 
has enabled a successful description of many important properties of strongly correlated 
systems in a unified frame. Such a success was accomplished since the DMFT fully 
includes the local quantum dynamics. On the other hand, the spatial correlations are 
treated in the mean-field level within the DMFT, and accordingly the DMFT is exact 
only in the limit of infinite dimensions. 
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Intensive efforts to improve the accuracy of the DMFT in finite dimensions have 
been made in many directions. Its cluster extensions [2] such as the cellular dynamical 
mean-field theory (CDMFT) [31 Hj, the dynamical cluster approximation |5J, and the 
variational cluster approach [6] include short-range correlations inside the cluster with 
the mean-field treatment of correlations between the clusters. The effects of spatial 
correlations are taken into account step by step as the cluster size is increased. The 
perturbative inclusion of the spatial correlations in a different direction has been tried 
by the diagrammatic expansions such as the dual-fermion method [7] and the dynamical 
vertex approximation [8j. The fluctuations due to spatial correlations are most severe 
in one dimension, leading to the expectation that one-dimensional interacting systems 
will provide strict tests for the efficiency of the attempts to incorporate the effects of 
spatial correlations. The Hubbard model is particularly useful since its exact solution 
is known in one dimension via the Bethe ansatz method [9]. 

The main purpose of our paper is twofold: We first examine systematically the 
improvement in the accuracy of the CDMFT with the increase of cluster and bath sizes. 
Secondly, we verify the performance of the CDMFT in describing the spectral properties 
of low- dimensional systems. For that purpose we investigate the properties of the one- 
dimensional Hubbard model at zero temperature within the CDMFT with an exact 
diagonalization method used as an impurity solver. The comparison of the CDMFT 
results with the exact ones obtained from Bethe ansatz solution shows that the CDMFT 
gives accurate results for the spectral gap and the energy density even with very small 
clusters used. Although for small local interactions odd clusters yield a metallic phase 
which is not present in one dimension, such a region is found to shrink quickly with 
the increase of the cluster size. Spectral weights computed from the CDMFT give a 
qualitatively good description of the spin-charge separation which is a characteristic of 
the one-dimensional interacting system. Furthermore they show excellent agreement 
with Bethe ansatz dispersions as well as the remarkable consistency in their intensity 
with the existing experimental and theoretical results, indicating the high accuracy of 
spectral properties within the CDMFT even with small clusters. Similar conclustions 
on the validity of the CDMFT in the one-dimensional half-filled Hubbard model were 
reached in a previous work [4] based on the investigation of spectral gaps. Here we 
emphasize that our analysis is extended significantly to the examination of the energy 
density and the spectral weights as well as the spectral gaps. Furthermore, the effects of 
the cluster size and the number of bath sites on the physical quantities are also clarified 
by the systematic study. 

This paper is organized as follows: section [2] gives a description of the one- 
dimensional Hubbard model and the method employed in this work. The results for 
the properties of the one-dimensional Hubbard model are presented in section [3j We 
give a summary in section HI 



Properties of the one- dimensional Hubbard model 



3 



2. Model and Method 

The Hubbard model is described by the Hamiltonian 

H = -t ( c L c ja + h- c + U Y n n n n ~ V Y 

(i,j}(7 i iu 

where c\ a (ci a ) creates (destroys) an electron with spin a at site i and the number operator 
is defined by n ia = c\ a Ci a . The hopping amplitude between nearest neighbors, the on-site 
Coulomb repulsion, and the chemical potential are denoted by t, U, and /i, respectively. 
We will consider the half-filled case /x = U/2 and express all the energy and the length 
scales in units of the hopping amplitude t and the lattice spacing a throughout the 
paper. 

We use the CDMFT to investigate the properties of the Hubbard model in one 
dimension. Within the CDMFT, the lattice problem is reduced to the interacting 
electrons on the cluster of N c sites hybridized with the noninteracting bath. The effective 
action S'eff for the cluster degrees of freedom is constructed by integrating out all the 
other degrees of freedom: 

f/3 



Jo 



fiver 



+ [ drf^Un^n^r), (2) 

where /i, v = 1, 2, • • • , N c are the indices of the sites inside the cluster and a Gaussian 
(dynamical) Weiss field Q~ l accounts for the effects of the other sites. The CDMFT 
imposes the self-consistency relation 

^{iun) - t c (iu n ) = Gy l(iu n ) (3) 

which determines the Gaussian dynamical Weiss field Q~ l . Here S c is the cluster self- 
energy which we obtain by solving the effective action in Eq. (T5]) and Gi oc is a local 
Green function for the cluster given by 

G Xoc (iu n ) = Y [(«^n + M)l - i(k) - t c (iuo n ) , (4) 
k 

where the momentum k runs over a reduced Brillouin zone of the superlattice composed 
of the clusters, t(k) is the Fourier transform of the hopping matrix, and the circumflex 
over symbols, " , represents an iV c x N c matrix form. 

In order to obtain the cluster self-energy we use a Hamiltonian formalism. We 
introduce an impurity Hamiltonian 

#imp = Yl E ^ C U C u<r + U Y1 n ^ n ^ 
fiver fi 

+ X^^IV + V X C U a la) + Yl e l^la a la, (5) 
fila la 
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where contains the hopping matrix elements inside the cluster and the chemical 
potential. The impurity Hamiltonian H imp contains additional degrees of freedom 
(I = 1, 2, • • • , Nb) on the bath which have energies q ct and are coupled with the cluster 
sites via V^. The Gaussian Weiss field corresponding to the impurity Hamiltonian in 
Eq. (ED is 

iuj n - ti a 



[Q l {iuJn; {tia, V^<t})W = ioJ n S^ - - . — vl ° , (6) 



I 

where 5^ u is the Kronecker delta. The parameters {e^, V^} are self-consistently 
determined by the relation in Eq. ([3]). 

The self-consistency loop starts from guessing the parameters {q, V^}. (The spin 
indices are omitted for a simpler notation.) We first compute the cluster Green function 
Gy,v by solving the impurity Hamiltonian through the exact diagonalization of an 
impurity Hamiltonian matrix. We then obtain the cluster self-energy from the Dyson's 
equation, 

t c = g- 1 - g-\ (7) 

which in turn yields the local Green function Gi oc in one dimension via 

n7v/Nc r ~ ~ ~ i- 1 dk 

{iuj n +fj)l-t{k)- S c (^ n ) —— . (8) 



G\ oc (iuj n ) 



TT/Nc 



2tt/N c 



The resulting local Green function G\ oc {iu n ) and the cluster self-energy Y, c (iuj n ) give a 
new Gaussian Weiss field Gnewi^n) 

Q-^{iu n ) = 6£{ia> n ) + S c (^ n ). (9) 

We extract new parameters {e;, V^} by minimizing the distance function y 2 

V max N c 

X 2 = -T Yl \lGnew( iuJ n))nv 

iV max -t- 1 n=Q ^ 

-[^(^{e,,^})]^! 2 . (10) 

The distance function is evaluated on Matsubara frequencies uj n = (2n + l)//3 with 
n = 0, 1, • • • , iV max , (3 = 200, and iV max = 800, and a conjugate gradient algorithm is 
used in the minimization. The self-consistency loop is repeated until the convergence is 
achieved. 

3. Results 

We first calculate the spectral gap S as a function of the interaction strength U to 
demonstrate the accuracy of the CDMFT. The gap is measured by the difference between 
the energy of the lowest electron peak and that of the highest hole peak in the single- 
particle density of states. We compare the results with an exact result given in a compact 
form 



U Ji smh(27ity/U) 
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Figure 1. Spectral gap 5 as a function of U for (a) even N c ; (b) odd N c . The gaps 
from the CDMFT calculation are denoted by symbols while the solid line represents an 
exact gap from the Bcthe ansatz solution. The empty symbols represent the spectral 
gaps obtained as U is decreased while those obtained with the increase of U are denoted 
by the filled symbols. The number of bath sites N\, is fixed to 6 for even N c and to 5 
for odd N c . 



which is obtained from the Bethe ansatz solution [TO]. The CDMFT results of the 
spectral gap in the one-dimensional Hubbard model are available in an earlier work pE]. 
In contrast to the previous work we present the results for various cluster sizes N c and 
numbers of bath sites Nf, which reveals clearly the roles of N c and Nb. 

Figure [1] shows the spectral gap for various cluster size N c . For the clusters with an 
even number N c of sites, the CDMFT yields a spectral gap very close to an exact value, 
as shown in Fig. [TJ(a). The two-site cluster calculation already exhibits a qualitatively 
correct behavior of a spectral g function of U, such as the appearance of a finite 

gap for every finite U. Indeed the deviation of N c = 2 results from the exact values is 
quite small; such deviation is further reduced by an increase of iV c . Unlike the earlier 
work |4j we have measured the spectral gap in the limit that the broadening of the peaks 
in the single-particle density of states vanishes. This limit discloses maximal numerical 
errors from the finite number of cluster or bath sites. We can see that the spectral gap 
within the CDMFT is in excellent agreement with an exact one in such a worst limit. 

Contrary to the case of even clusters, the agreement of spectral gaps from odd 
clusters (N c = 1, 3, 5) with exact results is not good, as shown in Fig. [JJ(b). They show 
a metallic phase over a finite region of small U before entering a Mott phase through 
a coexistence region, which is rather similar to the behavior in infinite dimensions. 
Although odd clusters exhibit poorer agreement with exact results than even clusters, 
the increase of N c also causes a rapid reduction in the critical interaction strength 
together with the shrink of a coexistence region. Indeed the cluster with N c = 3 
reproduces essentially an exact gap for U > 6. 

We also plot spectral gaps for two- and three-site clusters with various numbers of 
bath sites. As can be seen in Fig. [2J the spectral gap computed within the CDMFT is 
also dependent on the number of bath sites iV&. This arises because the Gaussian Weiss 
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Figure 2. Spectral gap S as a function of U for various Nf, with (a) N c = 2; (b) 
N c = 3. The empty symbols represent the spectral gaps obtained as U is decreased 
while those obtained with the increase of U are denoted by the filled symbols. The 
solid line represents an exact gap from the Bethe ansatz solution. 




Figure 3. Ground state energy density E/N as a function of U for (a) even clusters 
with Nb = 6, (b) odd clusters with Nb = 5 within the CDMFT calculation. The energy 
density with odd N c also shows a hysteresis in the coexistence region as the spectral 
gap does. Since the size of the hysteresis in the energy density is smaller than the 
symbol size in every system studied, we plot only the energy density obtained when 
U is increased. The exact energy density from the Bethe ansatz is represented by the 
solid line for comparison. 



field is approximated by a finite number of bath sites. Accordingly, the errors in the 
spectral gap decrease monotonically with N b , which is apparently faster for larger U. 

In addition to the spectral gap we also examine the ground state energy density in 
order to check whether such accurateness of the CDMFT shown in the study of spectral 
gaps persists in other physical quantities. The ground state energy density is calculated 



as 



E 

N 



^EE 







N r . 



-Tr 



H (k) + -S c 



x \{iuj n + i(k) - t c {iu n )]- 1 



(12) 
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Figure 4. Spectral gap i5 and ground state energy density E/N as a function of the 
cluster size iV c for various values of U. The results for U = 5,6, 8, 10 are denoted by 
the red squares, green circles, blue triangles, and pink diamonds, respectively. The 
lines represent exact values for the corresponding U with the same color. We used 
Nt = 5 for odd N c and Nb = 6 for even N c . 



where Hq includes kinetic and chemical potential terms in the Hamiltonian and the 
factor 2 comes from the spin degeneracy. The summation over Matsubara frequency 
is transformed to an integral in the zero-temperature limit. In Fig. [3] the energy 
densities within the CDMFT for various sizes of clusters are compared with an exact 
energy density. The exact energy density can also be calculated from the Bethe ansatz 
solution [9] given by 

E = T JpMJxMdu; _ u 

N J W [l + exp(iwtO] 4' 1 ] 

where J n (x) are the nth-order Bessel functions of the first kind. Except for the single-site 
calculation, all the calculations essentially reproduce the exact energy density. Figure @] 
displays the explicit dependence of the spectral gap and the ground state energy density 
on the cluster size N c . Both quantities show the even-odd oscillations with the increase 
of N c , and gradually approach the exact value in the whole range of U. 

The distance function \ 2 f° r even clusters is plotted on a logarithmic scale as a 
function of N b in Fig. [5j Since the infinite bath is projected into the space spanned 
by Nb sites in the exact diagonalization method, larger Nj, is expected to give smaller 
distance function, implying better convergence. In general the distance function \ 2 
decreases exponentially with increasing iV&. A small bump at Nb = 6 indicates that the 
baths of 4 and 8 sites are more favored than that of 6 sites in the calculation of the 
CDMFT. This is because it is difficult for 6-site bath to reflect the inversion symmetry 
and the particle-hole symmetry of the system. The effect of the cluster size N c on x 2 is 
negligible because the limitation of the accuracy in the distance function comes mainly 
from the use of a finite number of bath sites. 
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Figure 5. Distance function x 2 of the converged solutions for U = 
values of N c as a function of Nj,. 
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Figure 6. Local density of states p(ui) for various values of U for N c = 4 and TVf, 
We used a broadening factor e = 0.1. 



Within the CDMFT we next calculate the local density of states 

p(u) = --Im[TrG(w + is)}, (14) 

7T 

where e is a broadening factor. We plot the local density of states for various U in Fig. [6j 
The single band for U = is split to two bands, the upper and the lower Hubbard bands, 
with the local interaction U turned on. The two bands are located around ±(7/2. The 
Hubbard gap between the bands increases with the increase of U and becomes of the 
order of U in the strong interaction regime. It is interesting to note that the abrupt 
decrease of the density of states around the center of each Hubbard band separates two 
regions of the band, one with large density of state near the gap and the other with 
small density of states, which is more prominent for large U. The analysis of the spectral 
weights presented below provides an explanation for such a structure. 
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Figure 7. The colored plot of spectral weight A(k,w) for (a) U = 3.0, (b) U = 5.0 
and (c) U = 7.0. We used the cluster of size N c = 6 with the bath of Nt = 6 sites 
and computed A(k,u>) with a broadening factor e = 0.1. We normalized each weight 
by its maximal value. Some important dispersions and boundaries of the spinon-holon 
continuum, which are obtained from the Bethe ansatz solution, are also represented 
by solid lines, (see the text for details). 



The CDMFT allows us to compute the spectral weight A(k,u) as 

A(k, u) = —lmGi att (k, u + ie) (15) 

7T 

with the lattice Green function 

Here k is a vector in the original Brillouin zone and x M is the position of site [i in the 
cluster. The spectral weight shown in Fig. [7] gives an accurate description of the spin- 
charge separation which is a peculiar feature of the one- dimensional half-filled Hubbard 
model. When a hole or an electron is added to a one-dimensional interacting system, 
it is split into two collective modes, a spinon which is charge-neutral with spin 1/2 
and and a holon/antiholon which is spinless with a hole/electron charge. The Bethe 
ansatz method provides only the dispersions of such elementary excitations while their 
intensity could be calculated only in some special limits [TTI [T2] . Only recently has 



(16) 
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the reliable computation of the intensity been performed by dynamical density matrix 
renormalization group [T3l [H] and the cluster perturbation theory [15]. In Fig. [7] we 
marked some important branches: a spinon branch by A, a holon branch (0 < k < tt/2) 
and a secondary holon branch (n/2 < k < 7r) by B, the continuation of the holon branch 
from the negative-/c side by C, and the lower boundary of the spinon-holon continuum 
byD. 

The spectral weights computed for a 6-site cluster within the CDMFT exhibit an 
evolution with the variation of U which is quite consistent with Bethe ansatz dispersions. 
We can see clearly that spinon (A) and holon (B) branches show up at the locations 
predicted by the Bethe ansatz. Most weights are concentrated between the two branches, 
as has been observed in the experiments [TU [17] . As U is increased, we observe some 
weights between the branches B and C . Although some weights are found outside the 
branch C, they do not exceed the boundary imposed by the branch D. The fact that 
most spectral weights are located between the branches A and B explains the sudden 
jump in the density of states of the Hubbard bands which is mentioned above. For 
U = 7, for example, the large density of states in the region —4 < uj < —2 has the 
main contribution from the region between the spinon branch A and the holon branch 
B while the spectral weights outside the holon branch B form a region of small density 
of states up to the lower boundary of the spinon-holon continuum ( uj « —7 ). 

The plots of A(k, u) at some k points for U = 8 as a function of N c in Fig. [8] help 
us understand the distribution of spectral weights more clearly. The spectral weight 
at k = shows two prominent peaks, one around uj = —3 and the other around 
uj = —4.5, associated with spinon and holon branches, respectively. Other small peaks 
are found to form a continuum between the two peaks as N c is increased. At k = 7r/4 the 
system exhibits most spectral weight between the branches A and B with low-intensity 
continuum tails between the branches B and C. The spinon and holon branches merge 
at k = it/2, leaving one prominent peak around uj = —2 with a long tail up to the 
lower boundary of spinon-holon continuum. The spectral weights at k = 37r/4 and 7r 
also consistently show the continuum between the branches B and D. We find that 
the magnitudes of the spectral weights between the spinon and the holon branches are 
much larger than those in other regions. Our analysis verifies clearly that the CDMFT 
with a cluster of only 6 sites gives a quantitatively accurate description of the spectral 
properties in the one-dimensional half-filled Hubbard model. 

4. Summary 

We have investigated the properties of the one-dimensional Hubbard model in the half- 
filled case within the CDMFT. The exact diagonalization method has been chosen as 
an impurity solver. The spectral gap and the energy density which have been computed 
via the CDMFT even with very small clusters have been found to be in remarkable 
agreement with the exact values. The spin-charge separation, the inherent nature 
of the one- dimensional interacting systems, has also been observed in the CDMFT 
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Figure 8. Spectral weight A(k,u) for fc = 0, 7r/4, 7r/2, 37r/4, 7r and (N c ,N b ) = 
(2,8), (4,8), (6,6). We used a broadening factor £= 0.05. 



study of spectral weights. Furthermore, the CDMFT results for the spectral weights 
are in excellent agreement with Bethe ansatz dispersions and their intensity are quite 
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consistent with the existing experimental and theoretical results. It is not yet clear 
whether such remarkable accuracy in spectral weights can be achieved by the CDMFT 
with small clusters away from half filling, which would be an interesting future project. 
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